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We study the stability of anyonic models on lattices to perturbations. We establish 
a cluster expansion for the energy of the perturbed models and use it to study the 
stability of the models to local perturbations. We show that the spectral gap is 
stable when the model is defined on a sphere, so that there is no ground state 
degeneracy. We then consider the toric code Hamiltonian on a torus with a class of 
abelian perturbations and show that it is stable when the torus directions are taken 
to infinity simultaneously, and is unstable when the thin torus limit is taken. 

The presence of a gap is of crucial importance when considering the utility of topological 
phases for fault tolerant quantum computation [lj. A gap ensures exponential decay of 
correlations, the stability of operations to noise by localized perturbations and the possibility 
of adiabatically "dragging" anyons in order to perform gates. For a review of the ideas of 
topological quantum computing see e. g. [2]. Also, as shown by Hastings and Wen [3] the 
ground state degeneracy and structure of topological states (i.e. the exponentially small 
splitting of ground state degeneracy) is retained if a gap stays open along a path in the 
space of Hamiltonians connecting the perturbed and non perturbed Hamiltonian. Thus, it 
is important to verify the existence of a gap when searching for models with non abelian 
excitations which can be described by simple spin systems, and are rich enough to be useful 
for topological quantum computing. A recent example is Fendley's model [I], which is 
presumed to exhibit non-abelian quantum excitations, however, the existence of a gap is still 
an open issue. Understanding the presence and magnitude of gaps in generic Hamiltonians 
is also of critical importance for notion of adiabatic quantum computing [5]. 

In addition, temperature may play a competing role in destroying topological stability, 
and also affects the scaling of entanglement entropy as shown for the toric code and related 
Hamiltonians [HI El El E]- A way of addressing the fragility of the phase to temperature 
was proposed in [TU] and relies on confining the strings separating anyon excitations with a 
string tension terms. In [10] such terms are generated by coupling to a bosonic field. 

The presence of a gap in large quantum systems lies at the heart of our understanding 
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of correlations and phases. A common working assumption involved in treating many such 
problems is that the existence of a spectral gap ensures the stability of a phase to small 
perturbations. This is indeed the case whenever considering perturbations which are smaller 
in operator norm than the gap, by the convergence of perturbation theory. Thus a small 
localized perturbation cannot do real damage to the phase. However, such an argument is 
not useful when determining the stability of a given phase. Indeed, in any physical realization 
of a Hamiltonian describing a given phase, one may expect some degree of deviation from 
an idealized model, which is spread throughout the sample, and as such, may have a large 
norm, scaling as the volume of the system when a thermodynamic limit is taken. 

While showing that a system is gapless may sometimes be done by variational means as is, 
for example, done in the proof of Lieb-Schultz-Mattis theorem for spin 5 = 1/2 Heisenberg 
model [H] (and its generalization to higher dimensions p2]), establishing the presence of a 
gap may be substantially more involved, and cannot be addressed in this way satisfactory, 
as is seen from the closely related but much harder problem of establishing that a gap is 
present in the spin S = 1 Heisenberg chain (known as Haldane's conjecture). 

In view of the importance of a gap for the topological quantum computing idea, this paper 
is concerned with establishing a stability theory for such systems to quantum perturbations. 

Here we consider a class of Hamiltonians describing topological anyonic models on a 
lattice. An anyonic lattice model is a lattice A together with a finite set of anyonic charges 
C and fusion and braiding rules (which constitute a "unitary braided tensor category"). The 
Hilbert space T~C(A) of the theory is spanned by fusion diagrams. These describe the charges 
dj G C occupying sites i G A and the splitting process that was involved in creating them. 
We can write a state of the system in the form: 

0,1, Cl, X2, O2, C2, ..-XN-1, %-l,flJv) (1) 

Here one chooses some ordering of the lattice points, so that ai is the charge at site 1, 02 
the charge at site 2 and so on. These states may be conveniently represented by splitting 
diagrams as in Fig. [TJ The proportionality sign indicates that there is scalar factor relating 
the diagram and the normalized state. This factor is usually expressed in terms of the 
quantum dimension of the appearing charges, and is used to ensure the isotopy invariance 
of the various braiding and fusion operations that are carried out when expressed through 
these diagrams. The diagrams describe the generation of states from the vacuum. First, 
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FIG. 1: A basis element of the anyon lattice. 

creation of a pair a\,c\ = di in a channel Xi> then splitting c\ — > a 2 %2 through an allowed 
fusion channel with % 2 indexing the degeneracy. We continue this way until the charges at 
all points on the lattice are specified (note that some of the a^s may be vacuum, denoted 
a* = 1, corresponding to no charge at location i). We also assume that only a finite number 
of different charges are allowed to sit at a given site, and that these states are normalized. 

The fusion diagrams may be thought of as word-lines of particles which are created from 
the vacuum state. Closed diagrams correspond to processes which start in the vacuum 
state and return to it, and so correspond to a vacuum to vacuum amplitude, or a vacuum 
expectation value. The detailed algebraic structure of these theories was discovered and 
studied in [131 US HHJ US] ( see a l so [E3 HE] for a purely mathematical treatment). A useful 
summary of the essentials can be found in [HI [20], but a detailed treatment is not essential 
for the present work. Here the only important feature we use is that the vacuum expectation 
value of two braids which are unlinked is assumed to be factorizable: i.e. if the world lines 
consists of two closed braids C\ and c 2 which are unlinked (and for systems on torus, not 
completing a full cycle of the torus) then (c\ \)c<i) = (ci)(cz), where the expectation values 
are taken in a ground state of the system, i.e. a state containing no quasi particles. Thus 
the results of the present work may be extended to non-anyonic models once such a property 
can be established. 

Prime examples for anyonic lattices rise in the description of certain lattice spin Hamilto- 
nians such as Kitaev's toric code p] and the Levin- Wen models [21]. The latter supplies an 
explicit prescription for generating all discrete gauge theories and all doubled Chern-Simons 
theories in (2+1) dimensions. Another way of getting anyonic lattices is to consider at the 
outset an array of quasi particles with anyonic properties (such as the ones arising in frac- 
tional quantum hall effect), and assuming that their interaction can be completely specified 
by braiding and fusing operations. An explicit example for this is the array of Fibonacci 
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FIG. 2: A local perturbation Vi only allows braids and fusion in Ri x [0, 1], so that the end points 
sit on lattice cites. 

anyons known as the "golden chain" and its generalizations [221 123] ■ 

It is worth pointing out that even if the model is actually an effective description of a spin 
model, the particular lattice over which the corresponding anyonic lattice model is defined 
is not necessarily equivalent to the original lattice of spins, since particles may be associated 
with plaquettes of the original spin lattice and not with individual spins. 

As our unperturbed Hamiltonian, we take: 



hi\Xii °i) c i) X2, a 2, C2, ...Xn-i, a N-i, o,n) — h(ai)\xi, Oi, c 1; X2, a 2, ■ ■■Xn-i, aar-ij %) (3) 

where h is a function h : C — > R such that h(l) = and h(a,i) > h > 1 for any a, 7^ 1. It is 
clear from the definition that [hi, hj] = 0. We denote h m = maxc h(a). 

We now introduce the perturbations Vi to the system. Admissible V^s have a finite range 
and act in a compact simply connected region Ri C IR 2 , so that a Vi represents an action on 
a given state by combinations of braiding and fusing only charges aj on sites j G Ri, and 
such that the braid/fusion actions occur entirely within Ri x [0, 1] and the initial and final 
points of the picture lie on lattice points belonging to Ri (Fig. [2]). 

In general, even for neighboring sites Ri fl Rj 7^ 0, and the V^'s do not commute. 
However if Ri fl Rj — then [Vi, Vj] = 0, this can be verified directly from action on basis 
states (or for models arising from a spin lattice, by recalling that all local operators can be 
written as local products of spin operators), and illustrated in Fig. 13} 




(2) 



such that: 
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FIG. 3: V{S associated with disjoint R^s are commuting. 



Our first result is the following: 



Stability to quantum perturbations: 



Assume that the V^s are bounded: 



M W> < 1 e H{A), 



(4) 



the ground state of Hq is unique and the system is defined on a sphere (i.e. there are no non- 
trivial cycles). Then there exists a (3q > which is independent of the size of the lattice (but 
depends on the dimension and coordination number) such that for any (3 < @ , H +j3 J2ieA ^ 
has a unique ground state and spectral gap 7 > c(/3), where c(/3) is independent of lattice 
size. 

This part of the paper establishes the stability of the quantum phase described by H to 
quantum perturbations, when H has a unique ground state and may be viewed as a special 
case of [21] , where stability of classical models was studied using cluster expansions and the 
possible applicability to anyonic systems was taken into account, albeit without a concrete 
framework, since the importance of anyonic systems on lattices models and examples such 
as Kitaev's model and Levin- Wen models where absent. 

Next, we consider a situation of a system having a degenerate ground state, which is 
encountered in anyonic models, when placing the lattice on a manifold with a non-trivial 
genus. One of the simplest, and best known examples for a spin lattice Hamiltonian which 
yields an anyonic model is Kitaev's toric code pQ. This model was first proposed to be 
used as a quantum memory, by storing the quantum information in the subspace spanned 
by the ground states (the degeneracy of the ground states is 4 9 for a compact orientable 
surface of genus g). Here our motivation is to prove explicitly that Kitaev's argument for 
the stability of this memory does work. The argument is very simple: mixing between 
the ground states has to involve processes whereby a particle-antiparticle is created one of 
them tunnels around the torus and then they are annihilated, this cannot happen below Lth 
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order perturbation theory (where L is the length of shortest of the torus cycles), and so one 
expects an exponentially small splitting of the ground states o(e~ aL ) for some a. Much of 
what follows, is a way to establish this argument in a more explicit (and rigorous) way. 

Let us briefly review the toric code model. The system is composed of spins on a the 
edges of a square lattice. The Hamiltonian is given by : 

H TC = J2P+,s + J2 P ^P (5) 

s p 

where P +jS = |(1 — rijestarfs) a j )• Here star(s) is the set of four edges connected to a given 
vertex s on the lattice. P +)S projects on states where odd number of spins in a star(s) point 
in the x direction. Similarly P^ p — |(1 — rijgbouiidaryfp) a j) P r °j ec ts on states where an odd 
number of spins on the edges of a square plaquette p point in the z direction. One can easily 
verify that all terms in Htc mutually commute, so it can be diagonalized simultaneously 
together with all the projectors P +tS and Pn tP - The ground states can be thought of as states 
simultaneously satisfying the constraints P +jS = and Pu, P = for all stars and plaquettes. 

In this model there are two kinds of particles, denoted x, z. The x particles (termed 
"magnetic" in pQ) are introduced by applying 

x c = n < (6) 

where C is some set of edges of the square lattice to a ground state. Similarly z particles 
("electric") are introduced by applying 

to a ground state. In this system the z particles are associated with vertices touched by 
odd number of edges in C, while the x particles are associated with plaquettes surrounded 
by odd number of edges in C. Thus, for example given an edge a = an application of 
a z a to the ground state creates a pair of violations of the constraints = and P + j = 
and costs energy penalty 2. The set of charges contains now l,x,z, and any state in the 
system can be generated by successive generation of particles from the ground state, as in 
(OJ). Since the model is abelian, there is no degeneracy of the fusion channels and so the x 
labels are all trivial. 

When the model is placed on a sphere, the ground state is non-degenerate and so the 
stability theorem described above applies. We conclude that the toric code on the sphere is 
stable to perturbations. 
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Next, we consider the toric code on a torus, that is, our lattice of spins is now assumed 
to have periodic boundary condition, and we identify the sites = (i + Li,j) and 

(hj) — (hj + L2). The toric code has now four ground states, which may be characterized 
by applying loop operators to torus cycles. 

To understand what may happen under this degeneracy we proceed to study a particular 
kind of perturbations. The perturbations we consider are of the form: 



each Vi is composed solely out of sums and products of a z operators, and so can only create 
^-charges (electric). (Alternatively one may choose V% to be composed of o~ x operators and 
create "magnetic" charges). In j2S] a particular case of such a perturbation was shown 
numerically to preserve the topological phase. 
Our main results for this model are: 

1) One can completely destroy the gap if L 2 is kept fixed while L\ — > 00. 

2) If log(Li)/L 2 and log(L 2 )/Li are kept finite when the limit L 2 ,L 2 — > 00 then Htc is 
stable to the Vi perturbations. 

This is a preliminary result for a more general treatment of general anyonic systems on a 
torus (without the restriction on the type of perturbation). 



In this section we use techniques developed for the study of classical models and their 
perturbations and adapt them to the present setting. A classical model is a local Hamiltonian 
model defined on a tensor product Hilbert space (e.g. a spin lattice) in which the local 
terms are diagonalized simultaneously by a basis of product states. Indeed, while nontrivial 
topological Hamiltonians are never classical models, the particular form of the unperturbed 
Hamiltonians which we use, i.e. a commuting sum of local constraints, has much of the 
structure of a classical Hamiltonian. The stability of classical phases to temperature and 
quantum perturbations has been a rich field of study with well established results in the frame 
work of Pirogov-Sinai theory and it's subsequent refinements see e.g. [231 ESJ EH EH E2]. For 
classical Hamiltonians, stability of phases in the case of degenerate ground states usually 
depends on the Peierls condition. 




(8) 



I. 



PROOF OF STABILITY 
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In the present paper we consider the correction to ground state energy and gap of topo- 
logical lattice Hamiltonians by studying 



where \a) is a ground state of the system (possibly one of a set of ground states). If there 
is a gap, it follows that as iV — > oo 



where E a is the perturbed ground state energy in the cyclic subspace generated by applica- 
tions of if on |0), and 7 a is the gap above this ground state, in this cyclic space. 

It is convenient to compute log(ZAr) using cluster expansions. These have been extremely 
useful for the task of studying stability of classical models. However, the main difficulty in 
applying these considerations here, lies in the non-local nature of anyonic systems. This 
arises in two ways: 1) large clusters may be correlated with arbitrarily far other large 
clusters if they span non trivial loops around the torus, and as such cannot satisfy the usual 
factorizability properties needed to establish the convergence of the cluster expansion, and 
2) clusters which are disjoint but linked are not factorizable. In this section we consider 
systems defined on a sphere, so that correlation between arbitrarily far large clusters is not 
an issue and ground state degeneracy is not an issue. A system on a torus exhibiting a 
degenerate ground state and correlations between large clusters will be studied in the next 
section. 

There are a variety of ways of getting cluster expansions. Perhaps, the most straightfor- 
ward one is to write a Duhamel expansion of the partition function as is done in many works. 
However, here we chose to follow closely an elegant and more compact way by Yarotsky [30J , 
in order to get such an expansion, with the desired convergence properties. 

Given a set of vertices / C A, we associate to it a projection operator Pj = Ylsei 
which projects on the states in which no charge other then the vacuum 1 sits in I, and the 
operator Q = Yl seI P^~, which project on a state where all vertices in / carry a non-trivial 
charge. Using the inclusion-exclusion principle Z may be compactly written as: 




(9) 




(10) 




(11) 



c 



Here c = U^ =1 {Ik, Jk}, I C A, J C A\I, is a configuration on iV x 11, 
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N 



u(c) = (0\l[T Ih ,j h \0) (12) 



k=l 



and 



T ItJ = J2- llIHLle ~^ lhie ' HL Qj P ^ ( 13 ) 



Lcl 



Here I = \J ieI U J s t ieR D P ji i- e - ^he se ^ °f P°i n ts z's in A s.t. the site is affected by a 
perturbation in / and Hl = J2ieT^i + J2i^L ^i- 



To obtain the expansion (11) we have used that for any Jc A 



Yl Qj p ^\i\j = Id ( 14 ) 

JcA\/ 

where Id is the identity operator, and the inclusion exclusion principle through: 

e~ H = Y t ^2-l W ~ ]L] e-^ lhi e- Hl - (15) 

JCA Ld 

Therefore: 

e~»* = Y;UT IktJk (16) 

c fe=l 

The expectation value in u(c) is taken in the vacuum state. 

A configuration c can be identified with a set in A x (1, ..N). We define the support of a 
configuration by thickening the configuration in the time direction: define r{i, k) = (i,k + l), 
then supp(c) = c[J(rc) |J(r _1 c) in addition, we add a line for any two points in c which 
are next nearest neighbors on the lattice at a given "time slice", or are nearest neighbors 
in the "time direction" (Note that two different configurations can have the same support, 
depending on the interaction, for example if configuration (1) is the entire A x (1, ..N) and 
configuration (2) is A x (1, ..N) without one site they will have the same support the entire 
A x (1, ..N).). We denote the volume of a configuration |c| = of points in c}. 

The key property to evaluating is factorizability. For classical models, one usually 
shows that if c = c\ {Jc 2 , and supp(ci) fl supp(c 2 ) = then: uj(c) = u;(ci)u;(c 2 ). 

For anyonic models this is false. However we have a weaker kind of factorization property: 
If c = Ciljc2, and supp(c\) and suppfa) are non-linked, then: u(c) = u(cx)u)(c2). To 
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establish this property, we can take two configurations c\ and C2, and expand the T operators 



involved in (12) in terms of the perturbations Vj-s. The action of those on the ground state, 
in turn, can be expressed as a superposition of braid/fusion diagrams. It is clear that each 
term in such an expansion of c\ and C2 consists of a superposition braids, such that any pair 
of braids, one belonging to c\ and one to C2 are unlinked. By the axioms of the anyonic 
models we have that the vacuum expectation values of any such a pair factors. Thus we say 
that Ci and c 2 are compatible if their supports are not linked. 

Factorizability allows us to use the linked clusters in the following way. Write 



logZ iV = log Y[cu(c k ) = J2u(X) (17) 

irreducible&cunlinked ci,..c n X 

where X is a collection of polymers c& with multiplicities and, 

k ' Gx<G(X) 

Here G(X) is a graph obtained from X, which consists of n k vertices, where there are n k 
vertices which are associated to each c k in X . For two vertices associated with k% and &2, 
say, there is an edge if and only if c k \ and c k 2 are not compatible (i.e. their support intersects 
or links). The sum on the right is over complete connected subgraphs: G\ < G(X) i.e. G± 
is a subgraph of G(X), such that G\ is a connected and contains all vertices of G(X). In 
particular, if G(X) may be separated into two disconnected component, then there are no 
such G\, and so uj(X) = 0, which justifies the name "linked cluster expansion". 
To make this expansion useful we now need to establish two properties: 

1) Exponential bound on the number of clusters of given size and 

2) Exponential decay of u(c) with the volume of c. 



1. Exponential bound on the number of linked objects 

First, we prove that the number linked clusters of given volume L is bounded by v L for 
some v . This can be established in several ways but we use the following simple argument: 
Consider the doubled lattice (A x N) x {1, 2}. We can turn any linked configuration in ((Ax 
N), 1) into a unique connected configuration by adding segments of total length smaller then 
the linear size of the configuration (Fig. |1J) (since the distance between any two components 
which are linked to each other is bounded by their size). We then add this segments in the 
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((A x N),2) (one can think of adding it in "parallel") and connect them at some point to 
the original added segment as illustrated in Fig. |4j The volume of the added segments is 
then at most L + N where N is the number of connection points, therefore the size of the 
connected configuration in (A x N) x Z2 is bounded by 3L (see Fig. |4j). Since this mapping 
into the set of connected clusters in the doubled lattice is injective, we conclude that the 
number of linked clusters of size L is bounded by the number of connected clusters of size 
3L on the doubled lattice. The latter grows like v L for some constant v which depends on 
the dimension and connectivity of A [32] • 



2. Exponential bound on value of clusters 
Here we repeat the argument in [30]: Consider 

T hJ = J2-^ HLl ^^ lhl e- HUz ^ ZL) QjP A \r\j (19) 

Lcl 

where Hl(zi, ..Zl) = + SieL z iVi- Considered as a function of the Zj's, 

F t p(z 1 , ..z n ) = -ll'H^e-*^-*^) (20) 

Lcl 

is analytic for any ip. Note that if \zi\, \zz\, • ■ < 1, then for any tp, (ip\H L (zi, ..zl)\iP) < 
h m \I\ + f3\I\ < c\I\ for some c (independent of /3), where c > 1 is determined by the 
range of the V^s and h m . Therefore F(z±, ..z n ) < 2' 7 'e' 7 ' c and note that for any i E I, 
F<p(zi, z 2 , ..Zi-\, 0, Zi + x, z n ) = 0, since in this case the contribution from any set L C I\i to 



(20), is exactly cancelled by the contribution to (20) from L[Ji. By the Schwartz lemma 
we have that F^{z u ..z n ) < \ Zl \2^e^ c , for \z x \ < (3. It follows that F^{z u ..z n ) < 0%^°. 
In the same way, F^(zi, 0, Z3, ...) = 0, therefore again by Schwartz lemma F^,(zi,..z n ) < 
|^?2 |/52' 7 'e' 7 ' c , taking j^l < j3 and repeating this argument we find that: 

||^_ 1 |/|-W e -H i ^,.^)|| < p\i\ 2 \i\ e \i\° (21) 

Lcl 



We are left with evaluating We'^^^QjP^iW, but it is easy to see that each e~ hi act- 
ing on a state projected by Qj is bounded by e~ h °, and so we conclude that: HT^jH < 
(2/3e c )' 7 'e~' J ' h °. For any desired e picking (3 small enough and ho large enough will ensure 
that \u(c)\ < e\ supp ^\. 
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FIG. 4: For each linked configuration in the lattice A we can choose a corresponding connected 
configuration of larger size in Ai |J A2. No two configurations in A of length L correspond to the 
same configuration in Ai |J A2 

Having established these facts we can proceed by the usual means of cluster theory: The 
exponential bound on the weight of clusters together with the exponential bound on the 
number of clusters imply the absolute convergence of the linked cluster expansions if e is 
small enough. We identify the correction to the energy, and the gap by identifying the terms 



in the asymptotic expression (10). In the following, we sum over clusters which are defined 



on polymers on A x 00, note that clusters of time length l(X) > N cancel out when adding 



the terms below to get log Z N in ( 10 ) 



E 



X connected clusters starting at t=0 



(22) 



The o(e Nj ) term 



is: 



(KX) - N)ui(X) 

Xconnected clusters starting at t=0 with 1(X)>N 



(23) 



that this term is o(e~ Nl ) is a consequence of summing over clusters with volume larger then 
N, together with the exponential bounds on the number of clusters and their expectation 
values (The relation between 7 and e, v is described in [32])- Finally, the constant term is: 



C 



E KX)co(X) 

X connected clusters starting at t=0 



(24) 



In general, Eq, C depend on the lattice A and its volume. However, the exponent 7 is 
obtained irrespective of the size of A, therefore, the gap also holds when taking the thermo- 
dynamic limit. 



13 



To conclude we remark that the above treatment holds only in the cyclic subspace gener- 
ated by applications of H to the ground state. The treatment can be extended to the entire 
Hilbert space, by slightly perturbing the initial state and computing (0 + v \e~ NH \0 + v) and 
observing the changes to the cluster expansion. The resulting additional clusters are placed 
in the beginning or the end of uj(c) and it is clear that as such contribute as C\ + o(e~ jN ) 
to log Zn and do not ruin the gap [30J . 

II. TORIC CODE ON A TORUS 

We now turn to consider in more detail the toric code Hamiltonian Htc defined in Eq. 

The four ground states of Htc on a torus may be characterized by applying a z-loop 
operators to the torus cycles. These correspond to creating a particle antiparticle pair and 
then transporting the particle around the torus and annihilating it with the antiparticle. Let 
us denote by T\ and T 2 the two loop operators (one can choose, for example T\ = Ylfli °fi 
andT a = n£i*fi)- 

Here we deal with perturbations Vi which are functions of a z operators alone. Note that 
T\ and T 2 commute with Htc an d all of the Vi operators. The ground states can be chosen 
as eigenvectors of the T\ and T 2 . Since Tf — 1, their eigenvalues are ±1, we can build the 
projectors: |(1 ± T\) and =(1± T 2 ) which project on the symmetric and anti-symmetric 
subspaces. 

Let us denote \afj) where a, \i = ± the ground state of Htc such that Ti\afi) = <j\crfi) and 
T 2 \cr[i) = fi\<jfi). We stress again that any application of a z operators cannot mix between 
these states, therefore, the cyclic Hilbert space generated by applications of a z to these states 
is of the structure © CT =±,/Lt=±Wo^) where H afl = spcm{products of Vj applied to \a, /i)} 

A. Instability of thin tori 

An idealized description of Fractional Quantum Hall wave functions, which has been has 
been very successful is the "thin torus limit" [331 EH [35], [36] . In this situation one considers 
one of the directions of the torus infinitely long compared to the other. This yields a simple 
way of understanding the fusion rules and their application for various quantum hall states 
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by considering domain walls in a quasi Id object. 

Here we show that while the thin-torus limit is convenient from an algebraic/topological 
point of view, it also supplies an example of what can go wrong from the stability point of 
view: infinitesimal perturbations may ruin topological order if the ratio of the torus cycles 
is taken to infinity too fast. 

For the sake of exposition we look at the extreme situation where the torus has a small 
radius in one direction, so that we may consider a loop operator around this direction as 
"local". For thicker tori, loops wills arise at higher order of perturbation theory and one 
may show that an analogous situation to the one considered here happens. Thus, consider 
the perturbed toric code Hamiltonian: 

HTC + ^Y^lm (25) 

where l m = Yl^li °j, m - 

We note that [l m , H TC ] = 0, and that in the thermodynamic limit L 1 — > oo the pertur- 
bation is locally infinitesimally small. 

The ground state of this system is easily computed. Note that l m \crfi) = a\afi) so that due 
to the contribution from all the l m s the ground state subspace of Htc splits with energies 
— 1 and 1, and therefore there is no "topological" protection of the state, and it is not useful 
as a quantum memory anymore. It is worth noting that we can also make the system have 
a gapless like spectrum, by favoring a situation with a x-particle, or a domain wall present. 
This excitation can move around at a low energetic cost. Let us assume we have periodic 
boundary conditions of length 2Li, and take 

H = H TC -^-f](l m + l)-^- J2 ( 26 ) 

m=l m=Li+l 

Note that \a/j) are eigenstates of the system with eigenvalues —4. This system can be 
diagonalized exactly since all terms commute. If we introduce two x particles in the system 
in the form: \X,a/j) = Ylj^c a j \ a ^)^ where c is a cut from site 1 to site L 1 . We now 
note that: (1 + l m )\X, -/i) = 2\X, -fi) if m e [l,£i], while (1 - l m )\X, -/i) = 2\X, -fi) if 
m G [Li + 1, 2L\], since the string X can be commuted with the loop operator l m X = —Xl m 
for m E [l,Li]. Thus we have: H\X, —fi) = (2 — 8)|X, —fx) = — Q\X, —fx) this is the new 
ground state. The cost of moving the end of the string one step is very low: 4- if we move 
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in along the L\ direction (and if we move the string end in the L 2 direction). Therefore, 
as L 1 — > oo the spectrum becomes gapless. 

B. Torus 

As we have seen, the thin torus limit may be unstable to small perturbation. However, 
the situation is different if we take L\ and L 2 to the thermodynamic limit together. We 
find that if the ratio of L 2 and L\ doesn't blow up exponentially as we take that limit, then 
stability is retained. The idea we will use in order to show this property is to work in each 
of the four sectors (given the values of T 1? T 2 ) separately, and then compare the energy shifts 
and the gaps in each sector. 

We now compute the expressions for the energy and gap in each subspace as computed 
using the cluster expansion formalism from 

Z a ,(N) = (apt\e- NH M (27) 

To use the cluster expansion in this case we note that any cluster c containing a full 
loop around the torus, when expanded in the interactions Vj may contain loop operators. 
For general perturbations involving a z and a x operations, there will be no factorizability 
property of a pair of loops, even if they are arbitrarily far along the torus. This is clear from 
the following simple example: 

Li Li Li Li 

1 = HH^E"**^ + M n <iM M n <4M = 0, (28) 

i=l i=l i=l i=l 

where 1 < k < L\ is an arbitrary location for the second loop. The left hand side of this 
equation reflects that the loops fli^i a fi anc ^ Yli=i a ik are topologically equivalent. While 
the right hand side holds since 

M]Ja* k M = (aiM\T 2 f[a^T 2 \aiJL) = f] a^(-T 2 )T 2 M = ~M f] a^a/,). 

i=l i=l i=l i=l 

The crucial simplification when using only o z perturbations, is that any such loop will act 
on the ground states \afj) as a scalar and will commute with all other operators appearing 
in a cluster c, therefore factorizability is retained. We may now use the cluster expansion 
as derived in the previous section. 

However, we have to consider separately the subspaces with given numbers of x particles, 
since the Vi here cannot create those. Thus we consider first: 
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1. Energy shifts in the subspace with no x-particles present 

We proceed to compare the expressions for the energy obtained from the cluster expansion 
for each of the ground states. We write: 

Xconnected clusters starting at t=0 



The o(e Nj ) term 



is: 



^ (l(X) - N)u a „{X) (30) 

V connected clusters starting at t=0 with 1(X) >N 



where is given by (18) where expectation values are evaluated in the \a/j) state. 

We note that for any cluster of support less than a loop length, all ground states will have 

exactly the same contribution. Thus the difference between E ++ , E etc, is due to loops of 

support larger then mm(L 1; L 2 ). The contribution from such loops is of order exponential 
in the torus length i.e. the difference between say E + _ and E ++ is L 1 o(e _Ki2 ) for some k, 
where the Li factor comes from the number of possible starting points at t — 0. Since these 
states remained gapped from the rest of the spectrum by the same arguments as we had for 
the sphere, we conclude that the four new lowest energies are ensured to be exponentially 
close in energy if we demand that Lje~ KLj — > when Li,Lj — > oo. 



2. Energy when x-particles are present 

To complete the argument we must now consider the rest of the spectrum: Indeed, the Vj 
perturbations we introduced can only create one kind of particle, and so the considerations 
above only apply to the cyclic subspace containing no x particles. Here we have to show 
that the energetic penalty for introducing x particles, is larger then energy shifts in this 
subspace compared to the previous section (and contrary to the situation shown above for 
the thin torus). This can be addressed in the following way: we check the behavior of e~ NH 
on states containing a fixed number of x particles. 

To do so we introduce the x particles by applying to one of the ground states 

x c = n ( 3i ) 
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where C is some set of edges. Let us denote by n x the number of x particles associated with 
C. Since X c = 1 and X«t = X c it holds that: 

(X c e~ NH X c ) = ( e - NXcHXc ). (32) 



X c is a unitary transformation that acts as X c <J z a X c = —(J z a if a G C and X^cr^X^ = cr^ 
otherwise. Therefore XcHXc is a Hamiltonian where some of the a z operators have changed 
sign, and some of the Pq terms became Pq = 1 — Pq. These terms still commute with the 
rest of the Hamiltonian, and have eigenvalue 1 when applied to a state not containing x 
particles. The effect of Xc can therefore be summarized as: 

(Xce'^Xc) = ( e -*C*c+n.)) (33) 

where H c = H TC + Yl X c ViX c . 

We may compute the cluster expansion associated with the new Hamiltonian He- We 
have to identify the clusters which may feel the sign changes due to a z sign changes in some 
of the ViS. We note that the only clusters where a given a z appears an odd number of times 
(and so will be sensitive to the sign change) and which have non zero expectation in the 
ground state of H TC are those which contain closed loops (so they do not create particles). 
Such clusters are either encircling the entire torus, and so have weight o(e~ KL ), or have to 
go around one of the x particles. The contributions from clusters encircling a particular 
x charge can be bounded by YldLi c\de~ Kd < C2e~ K (This is equivalent to stating that the 
energy shift associated with a particular site is small when the perturbations are weak, or if 
the perturbation is translation invariant, the energy shift per unit volume is small). Thus 
if we demand that c 2 e~ K < e, we have that \E (H) — E (H C )\ < en x + o(e~ KL ), therefore, it 



follows from (33) that 

lo S ((Xce-»«X c »- log( n^ nAi _ £)+oi ^ L) (34) 

as X — > oo and we see that the states with x particles excited are separated by a gap of 
order n x (l — e) + o(e~ KL ) from the lowest states not containing x particles. 



III. SUMMARY 



In this paper we explored the role of perturbations on the behavior of certain topological 
lattice models. To simplify the treatment we started at the outset from the assumption 



18 



that the states are described by braid diagrams, with the property that any two non-linked 
braids factorize. We have proved that for weak enough perturbations, a topological phase 
defined by an anyon lattice on a sphere is stable. However, one may expect that more general 
systems may also be treated similarly. For systems on a torus we confined ourselves to the 
case of z perturbations of the toric code. Here we have shown how the stability of the model 
depends on the asymptotic ratio of the two torus directions. As expected, the system is 
stable whenever the ratio is not exponential. We expect such a result to hold for arbitrary 
local perturbations which are weak enough. This is the subject of an upcoming work. 
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